The Anderson-Mott transition as a random-field problem 
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The Anderson-Mott transition of disordered interacting electrons is shown to share many physical 
and technical features with classical random-field systems. A renormalization group study of an order 
parameter field theory for the Anderson-Mott transition shows that random-field terms appear at 
one-loop order. They lead to an upper critical dimension = 6 for this model. For d > 6 the 
critical behavior is mean-field like. For d < 6 an e-expansion yields exponents that coincide with 
those for the random-field Ising model. Implications of these results are discussed. 
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It is a well established result that electrons in a random 
potential at zero temperature undergo a metal-insulator 
transition as a function of the disorder, provided that the 
space dimension d > 2 . This is true for both noninter- 
acting 1^ and interacting [|j electrons, but the respective 
transitions are quite different in nature. For the Ander- 
son transition of noninteracting electrons no simple order 
parameter (OP) description, no upper critical dimension, 
and no Landau theory are known. The Anderson-Mott 
transition (AMT) of interacting electrons, on the other 
hand, has recently been shown to be conceptually sim- 
pler than the Anderson transition in that it allows for a 
simple OP description with the density of states (DOS) 
at the Fermi level as the OP. In this picture the metal 
(insulator), with a nonzero (vanishing) DOS at the Fermi 
level, represents the ordered (disordered) phase, respec- 
tively. This OP description leads to a finite upper critical 
dimension c?+, and to a Landau theory of the AMT 

This existence of an OP description raises important 
questions about the nature of the fluctuations that drive 
the AMT, and about its relation to other phase transi- 
tions in random systems. Consider a static, spin indepen- 
dent random potential u(x) that couples to the electron 
density. In a fermionic field theory this gives rise to a 
term in the action IHl 



dxM(x) Vn(x)-i/'n(x) , (1) 



with ip and "0 Grassmann fields, and n a Matsubara fre- 
quency index. Since the expectation value (Vi„(x)'(/;„(x)) 
determines the DOS, this means that the random poten- 
tial couples to the OP for the AMT. Magnetic transi- 
tion where a random field (RF) couples to the OP are 
known to have peculiar properties: The RF fluctuations 
are dominant over the thermal fluctuations, which leads 
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to d+ = 6 (rather than 4) |5| , and to a violation of hyper- 
scaUng even for d < 6 |^ . An e-expansion of the critical 
exponents about d = 6 leads to the famous 'dimensional 
reduction' problem and the critical dynamics have 
been proposed to be anomalous Q. It is then natural 
to ask whether similar phenomena are to be expected at 
the AMT. Physically it is plausible that interacting dis- 
ordered electron systems should display RF effects since 
they have the same type of frustration that occurs in RF 
magnets: The random potential favors a local electron 
arrangement that conflicts with the one favored by the 
electron-electron interaction. 

In this Letter we show that there is indeed a close 
analogy between the AMT and RF problems. Within a 
renormalization group (RG) treatment of the AMT we 
find that RF type terms appear which lead to d^ — 6. 
A 6 — e expansion to first order in e leads to critical 
exponents that are identical with those of the RF Ising 
model, and hyperscaling is violated due to a dangerous 
irrelevant variable. 

Our starting point is the nonlinear cr-model descrip- 
tion of interacting disordered electrons |9). This is a 
Gaussian field theory for a hermitian matrix field Q(x) 
with constraints ((5(x))^ = 1, with 1 the unit matrix, 
and trQ{x) = 0. Q is a classical field comprising two 
fermionic fields. It carries two Matsubara frequency in- 
dices n, m and two replica indices a, (3 (quenched disorder 
has been incorporated by means of the replica trick) . The 
matrix elements (5"m(^) general spin quaternions, 

with the quaternion degrees of freedom describing the 
particle-hole and particle-particle channel, respectively. 
For the sake of simplicity, in this paper we restrict our- 
selves to the particle-hole spin-singlet degrees of freedom, 
although the general model can be treated in the same 
way |lO|| . The matrix elements of Q can then be expanded 
as Q"l^ = I]r=o,3 rQuL^r with tq, 1,2,3 the quaternion 
basis. We write the action in the form 

+ 2H J dx tr (r!Q(x)) -^K,, [Q(x) o g(x)] . (2) 

Here G is a measure of the disorder, SI is a diagonal 
matrix whose elements are the Matsubara frequencies 
Un, and H is proportional to the free electron DOS. 
Ks < is a repulsive electron-electron interaction cou- 
pling constant in the particle-hole spin-singlet channel. 
We consider a short-ranged model interaction for simplic- 
ity. [Q o Q] denotes a product in frequency space which 
is given explicitly in Refs. HJ^. Notice that we have en- 
forced the constraint = 1 by means of an auxiliary, 
or ghost, matrix field A(x). 

The correlation functions of Q determine the physical 
quantities. Correlations of the Qnm with nm < deter- 
mine the diffusive modes which describe charge and spin 
diffusion, while the DOS is determined by (Q^^) It 
is therefore convenient to separate Q into blocks: 



dx tr 



A(x)[Q2(x)-l] + (a,,Q(x))' 
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Qlt = Q{nm)Q':±{x) + Q{n)Q{-m)q±{x) 

+e(-n)e(m)(g+)^f„(x). (3) 

We then integrate out the massless g-field. Since the 
action, Eq.(^), is quadratic in q this can be done exactly. 

We next expand Q and A about their expectation val- 
ues (rQl2(x)) = 5ro5i2Nm and (rAi2(x)) = 5ro5l2'ni> 

with 1 = (ni,ai), etc., 

rQl2(x) = 5rQ ^12 + r012(x) , (4a) 



rAl2(x) = 5rO 5l2 Ini + r'012(x) 



(4b) 



with (r0i2(x)) — (rV'i2(x)) = 0. Notice that Nn is pro- 
portional to the DOS at an energy a;„ measured from the 
Fermi surface |^ . The resulting action is quadratic in (j), 
but contains terms with arbitrary powers in -0. If one for- 
mally integrates out -0 one is left with an action in terms 
of the OP Q only. It can be shown that the terms of 
higher than second order in do not change the struc- 
ture of the resulting (/)-field theory It is therefore 
sufficient to integrate out ip in Gaussian approximation. 
This is easily accomplished with the result 



dx tr 



0(x)(-a^ + (A>)0(x) + g(((Q)0(x))N 
—u / (ixfr(/)*(x) 



(Q) 



2 J,2/ 



2u 

72d 



dx tr 



-— / dxtr 



AU^(x) + 



/2G 



(Q)0(x) 



dx (^htri (j){x 

i=>,< 

(Q)03(x)+0(x)(Q)02(x) 
2 



/2G 



dx tr 



Bcj){x) 



(5a) 



where we have scaled <^ by a factor of V2G. tr = J2ran denotes a trace over all discrete degrees of freedom, and 



htry^<c^ = X]rQSn>On<o dcnotcs 'half-traces' that extend only over positive and negative frequencies, respectively. 

[Q) and (A) , and are given by 

(5b) 



A and B are functions of (Q) and (A) , and are given by 

A((g),(A)) = (g)2-l + /((A)) , B{{Q),{K))^{K){Q)-2GH^l 
where f{{A)) is a matrix with elements r/12 = ^ro ^12 fni with 



fn 



G 



2'kTGK, 



n — m — l 



2ttTGKs 



(5c) 



for n > 0, and a similar expression for n < 0. The bare 
value of the coupling constant u is u — —G / {df / dl). The 
bare value of A is zero. However, a term of this structure 
is generated by the RG at one-loop order, and it is crucial 
to include it in the action. A comparison with Refs. 
shows that the half-trace terms have the characteristic 
replica structure of a RF term. Since it is quadratic in 
(j) this term contributes to the Gaussian (G) propagator. 
In the replica limit we find 



(G) 1 
r012(k)s034(p)) =(5(k + p)(5r 



1 



<5l3 (^24 + (-)'' <5l4 <^23 + 



4Ae(nin3) 



mi2 



16 /c2 

5rO ^12 <^34 



■mi2 



(6) 
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Here mi2 = [li+l2)l'i + u{Ni + N2Y. It is clear that the 
term proportional to A will increase the upper critical 
dimension by 2, as it does in the case of RF magnets. 
We therefore expect = 6 in the model given by Eqs. 

rather than d+ = 4 which one would conclude from a 
power-counting analysis of the action at zero-loop order, 
i.e. without the RF term 

Equation (|5^) requires some explanatory comments. 
(1) (Q) and (A) are determined by the conditions (</)) = 
and (V') = 0. At zero-loop order, these two conditions 
yield A{{Q), (A)) = B{{Q), (A)) = 0. This is the zero- 
loop order equation of state that has been discussed in 
Ref. jj]. It yields mean-field exponents, which constitute 
the exact critical behavior for c? > (i+ = 6. For d < 
6, renormalizations of the equation of state change the 
critical behavior. (2) In writing Eq. ( pa|) we have omitted 
some terms that are irrelevant by power counting for d > 
4. So are the terms of O (</>'') which we kept. However, as 
we will see, the latter couple to the RF coupling constant 
A and therefore must be kept. We have verified that none 
of the terms omitted, and no other terms generated by 
the RG, couple to A These considerations are in 

direct analogy to the case of a magnet in a RF ||] . 

We now perform a one-loop RG analysis of the action, 
Eq. , using a standard momentum-shell method . 
It is convenient to first consider the theory at criticality. 
Then we can put (Q) = (A) = 0, and consider the renor- 
malization of u for d = 6 — e. Since u is irrelevant for 
d > 4 we need to keep only contributions that are of order 
g = uA, A being relevant with a bare dimension of 2. A 
is not renormalized, and the d^-tcrm is not renormalized 
either to one-loop order, so the exponent ry = 0{e^). We 
obtain the following flow equation for g, 

sis + • 

with b the RG length rescaling factor. Eq. (|^) possesses 
a fixed point g* = 2e/9 + ©(e^). 

We now turn to the disordered phase, i.e. the insulator 
where N„=q vanishes and ln=o = ^ has a nonzero value. 
We thus put (Q) = 0, and renormalize the mass term (A) 
or I in the action. We obtain 

^^=2l-gl + Oig') , (8) 

and two equations that determine the renormalized equa- 
tion of state. 



where both N and / are to be considered as functions of 
fl. If we replace g in Eq. (^) by its fixed point value g*, 
we find that I scales like 
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We next invoke the equation of state to find the relation 
between I and the distance from the critical point, t. To 
this end, we expand the r.h.s. of Eq. ( |9a| ) for small values 
of I. The Z-independcnt contribution is t. At linear order 
in d = 6 one finds a term ~ and a term ~ ZlnZ. The 
prefactors of these two terms are related since df/dl — 
—G/u. Replacing g by g*, we can exponentiate and find. 

This holds for e > 0. For d > 6 one finds instead I ^ t a.s 
one would expect within mean-field theory. 

We can now combine Eqs. ( p^ ) and ( pT| ) to get the 
correlation length exponent v to first order in e. Since 
we also know 77 to this order, standard scaling arguments 
yield all other static exponents. We find 

- ^ + ^ + 0{e') , = + 0{e') , 7 = 1 + 1 + 0{e') + 0{e') ,6^3 + e + 0{e') . (12) 

In order to obtain /3 and S we have used the fact that 
u is dangerously irrelevant, so hyperscaling is violated 
(see , and the discussion below) , and have accordingly 
replaced d by d—2 in all d-dependent scaling laws. Notice 
that Eqs. (|l^) are identical with the corresponding results 
for a RF Ising model |,|l|l. 

We still need to determine the dynamical scaling ex- 
ponent z. For this purpose we obtain a relation between 



I, N, and from Eq. (9b). Expanding the r.h.s. for 



small I, going to criticality, and exponentiating, we find 
^^i+e/9+o(e^) ^ Q_ If combine this with Eq. (|ll|), in 
which we substitute t ^ N^/f\ then we find 

z = 3-e/2 + 0(e2) . (13) 

Notice that z — 5[3 / v = ijh, with yh the exponent of the 
field conjugate to the OP. This was to be expected since 
our RG did not involve any frequency integrals. 

Let us discuss our results. Grinstein has shown that 
hyperscaling is violated in RF magnets because because 
the quartic coupling constant is dangerously irrelevant. 
The same arguments apply here. In order to completely 
describe static scaling, we therefore need a third expo- 
nent, 9, in addition to the usual two independent expo- 
nents V and r/. 9 describes the fiow of u to zero, and 
enters all hyperscaling relations. Consider, for instance, 
the OP. Its scale dimension is d/2 —1-1- 77/2, which in 
the absence of dangerous irrelevant variables leads to the 
scaling law (3 — v{d — 2 + rj)/2. However, since u scales 
to zero like u ~ 6^^, and since N{u ^ 0) 1/-Ju, one 
has instead 

l3^v{d-2-e + 'q)/2 . (14) 

Our explicit one-loop calculation yields 9 = 2 -I- O(e^), 
but we expect 9 = 2 to all orders in e as is the case 
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in the RF Ising model |T^. Now consider the density 
susceptibiUty, dn/dfi, or the specific heat coefficient, 7 = 
limr-^o Cy/T, or the spin susceptibility, Xs- AH scale 
like an inverse volume times a time, so their naive scale 
dimension is d— z. The violation of hyperscaling changes 
this to d — 9 — z. Above we have seen that z is equal to 
Uh — Sp/v^ and therefore 

z={d-e + 2-r])/2 . (15a) 

Somewhat surprisingly for a quantum phase transition, 
the dynamics are not independent from the statics. This 
is due to the RF fluctuations being stronger than the 
quantum fluctuations. We thus find that all thermody- 
namic susceptibilities scale like the OP, viz. 

x{t,n)=b'^^+'^-'^-'^'^'^x[tb''\^b') , (15b) 

where x can stand for N , dn/dpL, 7, or Xs- ^ can stand 
for either external frequency, or temperature, or, in the 
case of TV, for the distance from the Fermi level. 

RF problems contain an anomalously divergent cor- 
relation function j^. In the present case this func- 
tion describes 'mesoscopic' fluctuations of the local DOS, 
C(x,y) = [iV(x)7V(y)]a,„ - N^, with N{yi) the (unaver- 
aged) local DOS at the Fermi level, and [..]av the ensem- 
ble average. At criticality C behaves like |g] C(fc ^ 0) ~ 
^-2+r/-e^ An experimental observation of such a strong 
divergence would indicate that RF features are indeed 
present at the AMT. 

We now consider transport properties. The charge or 
spin diffusivity D has a scale dimension of z — 2. Since 
no d is involved, 9 does not enter, and we have 

D{t,n) =b^-'' D{tb^/'',nb') . (16) 

One is also interested in the scaling behavior of the con- 
ductivity, (T = Ddn/dfji. The behavior of a depends on 
whether or not dn/d^ has an analytic background con- 
tribution in addition to the critical contribution given by 
Eq. ( |15b| ) . Ref. Q has argued that for the present model 
it does not. This yields 

ait,n) ^b^+''-'^a{tb^/'',nb'') . (17a) 

The conductivity exponent s is then 

s^i'{d-2-9) . (17b) 

In more general models dn/dfi might have a noncritical 
background contribution. In that case a will scale like 
the diffusivity, leading to 

s = iy{z^2)='^{d-2-9-i]) . (17c) 

In either case, Wegner scaling (i.e., s — v{d ~ 2) [p^), 
which in previous work in d = 2 -I- e had been found to 
hold for the AMT as well as for the Anderson transition 
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, is violated. This removes the requirement s > 2(d — 
2)/(i, which follows from Wegner scaling combined with 
the result of Ref. [|l6|, and has led to severe problems 
with the interpretation of certain experiments |^ . 

Finally, we note that our results imply that all of the 
complications (most of them not quite understood) that 
are known to occur in the RF magnet problem should be 
expected for the AMT as well. For instance, although 
9 ~ 2 to all orders in perturbation theory, this almost 
certainly changes due to nonperturbative effects 0| . d = 
4 is most likely some sort of a critical dimension. Finally, 
the non-power law dynamical critical behavior that has 
been proposed for RF magnets j|] should be expected to 
occur at the AMT as well, giving the AMT some aspects 
of a glass transition, despite the conventional power law 
scaling encountered in perturbation theory. 
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